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Abstract 

We provide an alternative proof of QIP=PSPACE to the recent breakthrough result |JJUW09[ . 
Unlike solving some semidefinite programs that captures the computational power of quan- 
tum interactive proofs, our method starts with one QlP-Complete problem which computes 
the diamond norm between two admissible quantum channels. The key observation is that we 
can convert the computation of the diamond norm into the computation of some equilibrium 
value. The later problem, different from the former semidefinite programs, is of better form, 
easier to solve and could be interesting for its own sake. The multiplicative weight update 
method is also applied to solve the equilibrium value problem, however, in a relatively sim- 
pler way than the one in the original proof IJJUW 09 ] . As a direct byproduct, we also provide 
a NC algorithm to compute the diamond norm of a class of quantum channels. Furthermore, 
we provide a generalized form of equilibrium value problems that can be solved in the same 
way as well as comparisons to semidefinite programs. 

1 Introduction 

The interactive proof system model, which extends the concept of efficient proof verification, has grad- 
ually become a fundamental notion in the theory of computational complexity since its introduc- 
tion MGMR85HBab85B in the mid 1980s. In this model, a computationally bounded verifier interacts 
with a prover with unbounded computational power in one or more rounds. The prover wants to 
convince the verifier to accept(reject) the input, and the verifier will make its decision based on the 
interacting process. 

The expressive power of this kind of interactive proof system model with at most polynomial 
rounds of communications is characterized [She92t ISha92t ILFKN92B by the well-known relation- 
ship 

IP=PSPACE 

through the technique commonly known as arithmetization. Many variants of the interactive proof 
system model have been studied by introducing new ingredients, such as the public-coin inter- 
active proofs IBab85tlBM88llGS89l . multi-prover interactive proofs HBOGKW88H , zero-knowledge 
interactive proofs fGMR85HGMW9lH as well as the competing-prover interactive proofs |1 K97|. 



*The work was completed when the author was visiting the Institute for Quantum Computing , University of Wa- 
terloo as a research assistant. 
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This paper mainly works with the quantum interactive proof system, which is defined IIKW001 
IWat99| in a similar way to ordinary interactive proof systems except the verifier and the prover 
have access to quantum computers. Similar to the classical cases, several variants of quantum in- 
teractive proof systems have been studied, including the ordinary quantum interactive proofs HWat99t 
IKWOOL public-coin quantum interactive proofs MMW05L zero-knowledge quantum interactive 
proofs HWat09allKob08llHKSZ08L multi-prover quantum interactive proofs BKKMV091IKM03H and 
the competing-prover quantum interactive proofs HGut05l IGW051 IGW07I1 . The complexity class 
QIP known as the problems having quantum interactive proof systems satisfies ilKWOOH 

PSPACE=IP C QIP C EXP 

Along with the introduction of the complexity class QIP, several complete problems for this com- 
plexity class have been discovered. The first complete problem, called close images, was first pro- 
posed in 2000 IIKWOOL Several relevant problems which can be reduced to close images were later 
discovered |RW05llRos08l . Especially, the quantum circuits distinguishability problem, which was 
proved MRW05II to be QlP-complete, serves as our start point to prove QIP=PSPACE. 

Recently, a big breakthrough |JJUW09| that proves QIP=PSPACE uses the primal-dual ap- 
proach 1IAK07H based on the multiplicative weights update method to solve a certain kind of semidef- 
inite programs that characterize the computational power of QMAM. The latter complexity class 
was proved [MVV05] to have equivalent expressive power as QIP. The multiplicative weights up- 
date method is a well-known framework (or meta-algorithm) which originates in many fields. Its 
matrix version, which was recently developed and discussed in a survey paper MAHK05all and 
the PhD thesis of Kale |Kal07l , was shown to be a great success in extending the potential appli- 
cations of this famous framework. Particularly, a combinatorial primal-dual approach for solving 
semidefinite programs (SDP) was proposed BAK07M based on the matrix multiplicative weights up- 
date (MMW) method. Under mild conditions, the primal-dual approach can be used to improve 
the time performance of many known approximation algorithms via semidefinite program relax- 
ations. The main advantage of the MMW method together with the primal-dual SDP solver for the 
purpose of simulating quantum complexity classes is that the resultant algorithm can be easily im- 
plemented efficiently in parallel^. By making use of the known result NC(poly)=PSPACE |Bor77l , 
we can solve these SDPs in PSPACE and hence show that some quantum complexity class is con- 
tained by PSPACE. Before the proof of QIP=PSPACE (JJUW09J, similar ideas were applied to show 
the containment QIP(2)C PSPACE 1JUW091 and QRG(1)C PSPACE |jW09). 



Unlike proving the result based on the formulation of the definition of the computational class 
QMAM, our proof starts with one QlP-Complete problem. The problem called quantum circuits 
distinguishability H computes the diamond norm between two mixed-state quantum circuits. If we 
were again to directly compute the diamond norm by its definition, we would encounter some 
SDPs or convex programs which are more complicated than the one we would have if start with 
QMAM. Although time-efficient algorithms have already been proposed I Wat09bl lBATS09ll to ap- 
proximate the diamond norm, it is unknown whether these methods can also be space-efficient, 
namely running in PSPACE. The crucial observation here, to circumvent the problem above, is to 
change the form of the diamond norm before the computation. The resultant problem (see The- 
orem Q} has a very neat form and can be expressed as an equilibrium value. Most importantly, 

1 We call any algorithm efficient in parallel if it is in NC. However, in our case, the size of the matrix representation 
of any quantum system is exponential in the input size. Thus, NC(poly) is considered as the final complexity class. 

2 It's been observed that one can also start with another QlP-Complete problem close images [KWOO] or the protocol 
to simulate QIP with competing provers in | GW05 1 to prove the same result in almost the same way. However, the 
particular choice here connects our algorithm to the computation of diamond norm which is of independent interest. 
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as we will see later, there is a space-efficient algorithm to solve the latter problem in PSAPCE. To 
our knowledge, this conversion for the first time establishes the connection between the computa- 
tion of the diamond norm and the computation of some equilibrium value. Precisely, Theorem [TJ 
claims that under certain conditions the gap between two promises of the diamond norm can be 
transferred to the gap between two promises of the constructed equilibrium value. 

The equilibrium value, or better known as the value with minimax-maximin form, is an impor- 
tant concept in theoretical computer science. For instance, many game theory related problems 
can be characterized naturally in this form. The fact that we can exchange the positions of min and 
max in any equilibrium value makes the problem well structured and provides a simpler iterative 
algorithm (also based on multiplicative weight update method) to approximate the equilibrium 
value than the one for SDPs. Similar ideas were known in the study of game theory^ before, and 
was applied [JW09| in the proof of QRG(1)CPSPACE. Due to the merits of the equilibrium value 
problem, the converted problem from the diamond norm in our paper has a relatively simpler 
solution (see Theorem H)) than the SDP considered previously [JJUW091 - As a result, our main 
theorem provides a simplified proof for the following fact. 

Corollary 1. QIP=PSPACE 

As a sequence of the connection we build between the computation of the diamond norm and 
the computation of some equilibrium value, we also demonstrate how our algorithm can be used 
to approximate the diamond norm to high precision of a class of channels efficiently in parallel. 
Precisely, we show 

Corollary 2. Given the classical description of any two admissible quantum channels Qq, Q\, the diamond 
norm of their difference \\Qq — Qi\\o can be approximated in NC with inverse poly-logarithm precision^ 

This result supplements the time-efficient algorithm for calculating diamond norm in MBATS091 
IWat09bL Although our algorithm only works for a special class of channels (also one of the most 
interesting cases), extensions of the current algorithm for a larger class of channels could be ob- 
tained if more complicated analysis is involved. 

It is interesting to compare the proof of QIP=PSPACE in this paper and the one in |JJUW09| . 
Our comparison represented here is threefold. First, obviously the two approaches diverge at 
the start point. However, this difference is actually subtle. As we mentioned before, the QIP- 
Complete we considered could be replaced by the close images problem almost with no change of 
the latter proof. If one investigates the result QIP=QMAM carefully, one will find this equivalence 
also comes from the close images problem. Recall that the diamond norm problem is more natu- 
rally formulated as some SDPs or convex programs just as QMAM does. It thus seems like we 
deliberately formulate the problem by an equilibrium value instead of a more natural formulation. 
We consider this as the main difference between the two approaches. 

Second, different formulations hence lead to the need of algorithms for different problems, 
SDPs and equilibrium value problems in our case. Due to the relation PSPACE=NC(poly), it suf- 
fices to find algorithms that are efficient in parallel. Fortunately, such algorithms for both problems 
can be obtained based on the matrix multiplicative weight update method. Nevertheless, the two 
algorithms are quite different in several aspects. We refer curious readers to Kale's thesis ||Kal07|| 
for complete details, while a brief comparison can be found below. We will refer the algorithm 

3 See the survey ]AHK05a| for more reference. 

4 Here the input size is exactly the size of those matrices representing two channels. Thus, the precision scales down 
to inverse poly-logarithm. 
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for SDP in Kale's thesis as the primal-dual SDP solver because there are indeed other methods 
for solving SDP also based on the matrix multiplicative weight update method. Conceptually, the 
primal-dual SDP solver exploits the duality between the primal and dual problems of a certain 
SDP while minimax-maximin equality is made use of for equilibrium value problems. It turns 
out the minimax relation gives a simpler proof of the correctness of the algorithm than the du- 
ality relation. Technically both algorithms require efficient implementation of some oracles. For 
those SDPs and equilibrium values about QIP, the oracle for the equilibrium value problem is eas- 
ier to design than the one for SDP| Furthermore, SDP solver faces an additional difficulty which 
is not applicable to the equilibrium value problem. As one restriction of the matrix multiplica- 
tive weight update method, any solution obtained for SDP problems only satisfies the constraints 
approximately. Namely, one needs to convert the raw solution into exactly satisfiable solution. 
However, there is no control in general about the change of the object function value after this 
conversion. Therefore, converting approximate solutions to exact satisfiable solutions without 
changing the object function value a lot is another difficulty in designing SDP solver. 

Finally, it is hard to compare the performance (e.g, in terms of time, space or other resources) of 
those methods for general equilibrium value problems and SDPs. We do not even know to what 
extent those methods can be applied to general equilibrium value problems or SDPs. The analysis 
might heavily depend on the particular form of the problem itself. However, some progress has 
been made recently BWulOtlGWlOB in finding efficient algorithms for a larger class of equilibrium 
value problems and SDPs. Particulary, there exists an equilibrium-value-based SDP solver HWulOl 
in addition to the primal-dual SDP solver. The new SDP solver provides a generic way to design 
efficient oracles, whereas a generic way of converting approximate solution to exactly satisfiable 
solution remains unknown. 

The rest of this paper is organized as follows. We briefly survey some preliminaries which 
will be useful in our proof in Section [2] The conversion from the QlP-Complete problem to some 
equilibrium value problem lies in Section|3l which is followed by the main proof of QIP=PSPACE 
in Section HI The algorithm for computing the diamond norm is discussed in Section [5] We con- 
clude the whole paper with the summary, Section where we provide further discussions about 
the equilibrium value problem and some open problems. Before the readers move on to the next 
section, there is one point to make clear. We will not take care of the precision issues with the NC 
implementation in the main part of this paper. Instead, we will assume such implementation can 
be made exactly and deal with precision issues in Appendix [B] 

2 Preliminaries 

This section contains a summary of the fundamental notations about the useful linear-algebra facts 
in quantum information. For the most part of this section, it is meant to make clear the notations 
and the terminology used in this paper. For those readers who are not familiar with these concepts, 
we recommend them to refer to IBna97l IKSV021 INC001 IWatOHH . 

A quantum register refers to a collection of qubits, usually represented by a complex Euclidean 
spaces of the form X = C E where E refers to some finite non-empty set of the possible states. 

5 It is not easy to compare directly since the algorithm in |JJUW09| unpacks everything and only uses one sub- 
routine, namely projection onto positive subspace (the same as the one in our algorithm). Nevertheless, if one rewrites the 
algorithm component by component, one could find out that oracle is slightly harder to solve. Moreover, additional 
assumptions like the invertibility of some matrices are also necessary to solve that oracle. 
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For any two complex Euclidean spaces X ,y, let L (X ,y) denote the space of all linear map- 
pings (or operators) from X to y (L (X) short for L (X, X)). An operator A G L (X, y) is a linear 
isometry if A* A = \x where A* denotes the adjoint (or conjugate transpose) of A. 

An operator A G L (X) is Hermitian if A = A*. The eigenvalues of a Hermitian operator are 
always real. For n = dim X , we write 

A a (A) > X 2 (A) > ••• > A„(A) 

to denote the eigenvalues of A sorted from largest to smallest. An operator P G L (X) is positive 
semidefinite, the set of which is denoted by Pos (A"), if P is Hermitian and all of its eigenvalues are 
nonnegative, namely A„(P) > 0. An operator p G Pos (X) is a density operator, the set of which is 
denoted by D (X) , if it has trace equal to 1. It should be noticed that a quantum state of a quantum 
register X is represented by a density operator p G D (X). 
The Hilbert-Schmidt inner product on L (X) is defined by 

(A, B) — TrA*B 

for all A, B G L (X). 

A super-operator (or quantum channel) is a linear mapping of the form 

Y : L (X) — > L (y) 

A super-operator Y is said to be positive if Y(X) G Pos (3^) for any choice of X G Pos (X), and is 
completely positive if Y ® H-l(z) is positive for any choice of a complex vector space Z. The super- 
operator Y is said to be trace-preserving if Tr Y(X) = Tr X for all X G L (A"). A super-operator Y is 
admissible if it is completely positive and trace-preserving. Admissible super-operators represent 
the discrete-time changes in quantum systems that, in principle, can be physically realized. 
One can also define the adjoint super-operator of Y, denoted by 

Y* :L(y)^L{X) 

to be the unique linear mapping that satisfies, 

<B,Y(A)) = (Y*(B),A) 

for all operators A G L (X) and B G L (y). 

The Stinespring representations of super-operators is as follows. For any super-operator Y, there 
is some auxiliary space Z and A, B G L (X ,y <g) Z) such that 

Y(X) = Tr z AXB* 

for all X G L (X). When Y is admissiable, we have A = B and A is a linear isometry. 

A quantum circuit is an acyclic network of quantum gates connected by wires. The quantum 
gates represent feasible quantum operations, involving constant numbers of qubits. In a mixed 
state quantum circuit MAKN98L instead of using unitary operations as quantum gates, we allow the 
gates to be from any set of quantum admissible operations. In this more flexible circuit model, 
some part of the qubits might be discarded (or traced out) during the evolution of the circuit. 

The trace norm of an operator A G L (X) is denoted by || A \\\ and defined to be 

|| A ||i = Tr VA*A 
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When A is Hermitian, we have 

|| =max{(P -Pi,A) : P 0r P 1 e Pos (X) ,P + P x = l x } 
The diamond norm of a super-operator Y : L (A?) — > L(y) is defined to be 

||Y|| = max ||Y®l*(X)||i 

l|X||i<l 

Because of including the effect of using entanglement between the input and some auxiliary space, 
the diamond norm serves as a good measure of the distinguishability between quantum opera- 
tions. Furthermore, we can show 

Fact 1. MKSV02H If a quantum channel Y can be represented by Y(X) = Tr^ AXB* where A, B £ 
L (X, y ® Z), define the new channels 

Y A (X) = Try AX A* , Y B (X) = Try BXB* 

then the diamond norm of this channel Y is equal to 

||Y||o =F ma x(Y A ,Y B ) (1) 

where 

F m ax(Y A ,Y B ) =max{F(Y A ( (? ),Y B (0) : Q,£eD(X)} 

and 

Hp,Q) = \\VpVq\\i 

which is a generalization of the fidelity between quantum states. 

3 Conversion of the QCD Problems to Equilibrium Value Problems 

Definition 1 (Quantum Circuit Distinguishability). For any constant a, b, such that < b < a < 2. 
We define a promise problem QCD„ b as follows. Given the description of any two mixed-state 
quantum circuits Qo and Q\, which are admissible quantum channels from L {X ) to L (y), exactly 
one of the following conditions will hold: 

1. \\Qo-Qi\U>a 

2. HQo-Qillo < b 

QCD fl y will accept on the condition (1) and reject otherwise. 

It was proved by Rosgen et at MRW05B that for any constant < e < 1, QCD 2 _ £ /E is QIP- 
Complete. A careful reformulation of this problem will enable us to rephrase this promise problem 
in term of an equilibrium value problem. 

Assume there exists some space Z ® Q that will be constructed later, we define a min-max 
value A(S) to be min p <=D(X ®x 1 ) max ner (n,H(|0)) where S is a linear super operator mapping 
from L (X <g> X x ) to L (Z <g> Q) and T = {n : < n < 1 z ®q}- The x o, X\ in the above definition 
are isomorphic copies of X . Further investigation shows the value A(H) is also an equilibrium value. 
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Given D (Xo ® X\) and T are convex and compact sets and the (IT, E(p)) is a bilinear func- 
tion over them, it follows from the well-known extensions of von' Neumann's Min-Max Theorem 
llvN28llFan53ll that 

A(E) = min max (TI, E(p)) = max min (IT,E(p)) (2) 

peDiXg^X^ ner r// ner p<ED(Xa®X x ) 

The equilibrium value A(E) is the quantity represented by the two sides of the above equation. 
Furthermore, any pair (p, Yl) which makes the function reach the equilibrium value is called the 
equilibrium point; or, equivalently, that 

min (fLS(p)} = <fl,S(p)> = max (n,E(p)) 

The existence of the equilibrium point follows easily from Equation 121. Careful readers might 
notice the equilibrium value's form defined in this paper is related but slightly different from the 
one defined in the proof of QRG(l)cPSPACE | [JW09| . In the latter's definition I is the set of the 
density operators. Thus, the equilibrium value will be the largest eigenvalue or Ceo norm in some 
sense. However, our definition of T makes the equilibrium value be the summation of all positive 
eigenvalues. Moreover, in the situation of later discussion in this paper, the equilibrium value 
turns to be half the C\ norm. 

Our main theorem of this part says the two promises of any QCD problem, or equivalently 
of any diamond norm of the difference of two admissible channels, can be transferred to the two 
promises of the value of A(E) where E can be constructed efficiently from the input to that QCD 
problem. Precisely, 

Theorem 1. For any instance of the QCD a b problem, there exist some space Z ® Q and a linear super 
operator Efrom L (Xq ® X\) to~L(Z ® Q) where the space Xq, X\ are isomorphic copies of the space X 
such that 

f X(S) < ||Q -Qi||o>«; 

[ A(E) > ^, ||Qo - Qi||o < 

where A(H) is the equilibrium value defined above. Further more, such a super operator S can be constructed 
efficiently in parallel from the input to the QCD n /b problem. 

Before we get into the proof of the theorem, it might be helpful to see where this theorem leads 
us to. Since for any < e < 1 the QCD 2 _ £/E is QlP-Complete, we can choose a constant e' such 

that A(E) is either at least or at most ^^'- £ ' 2 where — s/ ' 4:£ ' 2 ^ £ ' 2 and there is a constant gap 
between the two promises. For example, if we choose e = 0.1, then the two promises become 

either A(E) > 0.95 (namely, ||<D||o < 0.1) or A(E) < 0.32 (namely, ||<D||o > 1.9) 

Thus, in order to simulate QIP, it suffices to compute A(S) approximately to distinguish between 
the two promises. This accomplishes the conversion we need for the next step of the whole proof. 
A simple proof for Theorem[I]is available below. 

Proof. For any instance of QCD a b , we are given the classical descriptions of two mixed-state quan- 
tum circuits Qo and Q\, which are admissable quantum channels from L [X) to L (y). Thus, we 
could describe the two circuits using the Stinespring representation of quantum channels. Pre- 
cisely, 

Qo(X) = Jr z (A XA* ), Q X (X) = Tr Z {A X XA\) 
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where Z is the auxiliary space and Aq, A\ G L (X, y ® Z) are linear isometries. Then, we have, 

O(X) = Q (X) - Qi(X) = Tr 2 (A XAS - ^XAJ) = Tr zm (2C Q XC\) 
where Q is a complex Euclidean space of dimension 2 and Q, Q G L (A", 3^ <g> Z ® Q). Moreover, 

^Mt)' <?=^(*-*) 

It is easy to see that CgQ = Q = 1^- given that AqAq = A\A\ = 1%. To compute the diamond 
norm of O, we define 

O a (X) = Tr^QXq), <D B (X) = Tr^QXC*) 

Due to Fact[T], we have 

|| O ||o =F max (20 A ,20 B ) (3) 

It is interesting and useful to see that we can use one density operator pGD (Xq ® X\) to represent 
q, Z, G D (X) where the space Xq, X\ are isomorphic copies of the space X. 
Precisely, let 

& A (X) = <S> A (Tr Xl (X)), $b(X) = <D B (Tr* (X)) 

It is obvious that 0^,<3> B are admissible quantum channels from L (Xq ® X\) to L [Z ® Q). 
Define 

F~ (20^,20^) = max{F(2<D A (p) 7 2<D fl ( ( o)) :p£D(X ® X t )} (4) 

By taking p = q ® £ in the Equation ffl, we have F^ x (20^,20^) > F max (20 A ,20 B ) . To 
see the reverse side, we can take q = Tr^p and £ = Tr^ oj o in Equation (3). Thus, we have 
F~ (20^,20^) = F max (20 A ,20 B ). Namely, 

||<E>||o = max{F(20l(S,20 B V)) :p£V(X ® X{)} 

Since in QCD problem we have the promise that either 1 1 <3E> 1 1 > # or ||<&||o < b. Due to the 
Fuchs-van de Graaf Inequalities, for any G D (X) , 



i-\\w-ai<HQ,o<\j^-\\w-a\ (5) 

and let T = {n : < n < 1 z ®q}. By substituting ¥(® A , $b) into Eq and make use of the fact 
F(20 A ,20 B ) = 2F(<3> a ,O b ), then we have when ||<3>||o > a, 

min max ( n, <I>a (p) -<E>b(|0)> = min - B (p) ||i < 



peD^o®^) ner \ / p<eT>{Xk®x x ) 2 2 
and when || O || < b, 

min max/n,^^) = min - <^(p) ||i > 

l oeD(A'o®A' 1 ) ner \ / peD^o^^) 2 2 
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1. Initialization: Pick a fixed e < \, and let =l^eL (Af), N = dim A". 

2. Repeat for each t = 1, . . . ,T: 

(a) Let the density operator = / Tr WW 

(b) Observe the loss matrix G L (A") which satisfies -1* < < or < <l x , 
update the weight matrix as follows: 

W( f+1 ) =exp{-sf^M { ^) 

T=l 



Figure 1: The Matrix Multiplicative Weights Update method. 

Let H = <$>a — Og and A(H) be the equilibrium value defined before. Finally, we have 

A(S)<vCT, ||Qo-Qi||o>fl; 
A(E)>^, ||Q -Qi||o<b. 

As we can see through the proof, the desired super operator S is constructed explicitly from the 
input circuits Qo, Q\. Moreover, every step in the construction only involves fundamental oper- 
ation of matrices. Due to the facts in Section 14.21 we are able to construct such S efficiently in 
parallel. □ 

4 Multiplicative Weights Update method for Computing Equilibrium 
Values 

The multiplicative weights update method introduced in Section[l]is a framework for algorithm design 
(or meta-algorithm) that works as the one shown in Fig [TJ This kind of framework involves lots 
of technical details and we refer the curious reader to the survey and the PhD thesis mentioned in 
the introduction. However, for the sake of completeness, we provide the main result which will 
be useful in our proof. It should be noticed that {M^} is the freedom we have in this framework. 

Theorem 2. After T rounds, the algorithm in Fig\l\guarantees that, for any p* G D {X), we have 

(1 - e) £ M») + (1 + e) £ U\ M«) < lp\ £ M« \ + ^ (6) 

>o <o \ t=i / 6 

Here, the subscripts > or < in the summation refer to the rounds t where < M( f ) < 1 or —1 < 
Ifiyl < o respectively. 

Since in our consideration, it always holds that < M( f ' < 1. It suffices for our purpose to 
prove a simpler version of the theorem |2] although the proof is almost the same as the one for the 
general version. 
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1. Lete = I and T 



161nN 

S 2 



. Also let = l x , N = dim (X). 



2. Repeat for each t = 1, . . . , T: 

(a) Let = W (t ) / Tr and compute the S(pW ). Then let be the projection onto the 
positive eigenspaces of S(pW). 

(b) LetM (t) = (E*(IlW) +l x )/2, and update the weightmatrix as follows: 



W( m > = exp(-ef^M^) 

T=l 

3. Return i X^Li ( FI^, E(pW ) ) as the approximation of A(E) . 



Figure 2: An algorithm that computes the approximation A(S) with precision 5. 

Theorem 3. Assume < My' < 11 for all t, after T rounds, the algorithm in Fig\l\guarantees that, for 
any p* G D {X), we have 

(l-e)t (p {t \M^\ < (V, t mA + ^ (7) 

We put off the proof in the appendix part. It will be sufficient to just remember this theorem in 
the first reading and skip the details. 

4.1 Algorithm for Computing Equilibrium Values 

Using the multiplicative weight update method to compute some kind of equilibrium values was 
known before, for instance the equilibrium value of zero-sum game (an algorithm to compute this 
value can be found in Kale's thesis ||Kal07|| and the survey ||AHK05bL The reference for similar 
algorithms with different purposes can be found in the survey ||AHK05bll ). However, to compute 
the equilibrium value defined in our form, we need to adapt the old idea to the new situation. 

In order to compute the equilibrium value A(H), we design an algorithm as shown in Fig|2j 
This algorithm takes the descriptions of the two mixed-state quantum circuits as input, and then 
compute the S = — in Theorem [TJ and output the equilibrium value A(E) with precision 5. 
Namely, the return value A satisfies | A — A(H) | < 5. 

Before proving the correctness of the algorithm, one might want to compare the algorithms in 
both Fig [T] and Fig [2] The only change in our algorithm is that we propose a way of computing 
for each round t. As we mentioned before, {M^ } is the freedom we have in this framework. 
Different designs of {M"} can lead to different applications of this framework. For instance, 
the primal-dual approach for semidefinite programs in Kale's thesis [Kal07] is an example of the 
design of {M^} that provides a good application. 

Theorem 4. Using T = 16 ^ N rounds, the algorithm in Fig^returns the approximated value ofX(E) 
with precision 5. Namely, we have the return value A satisfying 

|A-A(S)| < 5 
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Proof. First note that for any TL^> computed during the process, 

Vp€D(#), |/p,S*(n«)\| <1 



since E(p) is the difference between two density operators. Thus, MW = (E*(n( f ') +l^)/2 
satisfies < < 1*. 

Then apply Theorem [H we have, 

(l-e)E (pW,MW) < lp\ E MM) + — (8) 

T=l \ T=l / ^ 

for any density operator p* G D (A"). Substitute = (E*(nM) + t x )/2 into Eq J3 and divide 
both side by T, note that (p^\ ) < 1, then we have 



I E ( P w # s*(nW)) < 1 U, E E ( (nW) +2, + ^ 



By choosing £ = | and T 



161nN 



, we have 



a = I E (pW a*(nW)) < 1 ( P *, E s*(n< T >) ) + * 



1 T=l \ T=l 



(9) 



In each step, IlW is returned as the solution to maximize y n^, S(p W ) ^ for any fixed pW . Due 
to the definition of the equilibrium value in Eq |2|, the equilibrium value A(S) < (Tl^\ E(p( f ')^) 
for any returned n( f '. On the other side, choose (p, ft) to be any equilibrium point and let p* = p, 
then we have (p,S*(nM)\ < A(E). Using inequality 0, we have 



A(S) < A < I ( p, E S*(n«) ) + <5 < A(S) + 5 (10) 



r 

E 

± \ t=i / 

which completes the proof □ 



To distinguish between the two promises in Theorem [IJ we let 5 = 0.2 and make use of the 
approximated equilibrium value returned in the algorithm. If the value is closer to 0.95, then it is 
the case that || <E> ||<> < 0.1. Otherwise, it is the case that H^Ho > 1.9. Thus, we solve the promise 
QCD problem in this way. 

4.2 Simulation by bounded-depth Boolean circuits 

We denote by NC the class of promise problems computed by the logarithmic-space uniform 
Boolean circuits with poly-logarithmic depth. Furthermore, we denote by NC(poly) the class of 
promise problems computed by the polynomial-space uniform Boolean circuits with polynomial 
depth. Since it holds that NC(poly)=PSPACE, thus in order to simulate the algorithm above in 
PSPACE, it suffices to prove that we can simulate the algorithm in NC(poly). 
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There are a few facts about these classes which are useful in our discussion. The first fact is 
the functions in these classes compose nicely. It is clear that if / G NC(poly) and g G NC, then 
their composition g o f is in NC(poly), which follows from the most obvious way of composing 
the families of circuits. Another useful fact is that many computations involving matrices can be 
performed by NC algorithms (Please refer to the survey HGat93H which describes NC algorithms 
for these tasks). Especially we will make use of the fact that matrix exponentials and positive 
eigenspace projections can be approximated to high precision in NC. A more careful treatment on 
those issues can be found in Appendix [Bj 

Since we are able to perform matrix operations with sufficient accuracy in NC, it remains to 
show the ability to compose all the operations in NC(poly) and thus in PSPACE. Precisely 

Theorem 5. The algorithm shown above can solve QCD problems in NC(poly), and thus in PSPACE. 

Proof. To simulate the algorithm, it suffices to compose the following families of Boolean circuits. 

1. A family of Boolean circuits that output the representation of the quantum channel S (in 
Theorem [T]) generated from the input x, namely, the descriptions of two mixed quantum 
circuits. 

2. Follow the algorithm in Figure |2l Compose all the operations in each iteration. Consider 
the fact that fundamental matrix operations can be done in NC and the number of iterations 

is polynomial in the size of x since 5 is a constant and N is exponential in the 



T 



161nN 



S 1 

size of x. 

3. The circuits to distinguish between the two promises by making use of the value returned in 
the circuits above. 

The first family is easily done in NC(poly), by computing the product of a polynomial number of 
exponential-size matrices which corresponds to the mixed quantum circuits. The second family 
is in NC by composing polynomial number of NC circuits. The third one is obviously in NC. 
The whole process is in NC(poly) by composing the NC(poly) and NC circuits above, and thus in 
PSPACE. □ 

It follows from the Theorem [5|that we can solve QCD\ 9 0.1 problems in PSPACE. Since QCDi 9 0.1 
is QlP-Complete problem, and any polynomial reduction to that problem can be easily done in 
NC(poly) by computing the product of a polynomial number of exponential-size matrices, thus 
we have QIPCPSPACE. Combining with the known result PSPACE =IPCQIP, we have, 

Corollary 1. QIP=PSPACE 

We notice that all the proof above so far is based on the assumption that all the matrix oper- 
ations can be simulated exactly. However in practice, we will need to truncate the precision to 
some place for some operations to be performed. Fortunately, this won't be an essential obstacle 
for the implementation of the algorithm. As mentioned in the introduction, all those issues will 
be handled in Appendix|B]without any change of the main result. 



5 Algorithm for computing the diamond norm 

Now it is our turn to discuss the computation of the diamond norm for a special class of quantum 
channels. Consider the following promised version of the problem first. 
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Definition 2 (Promised Diamond Norm Problem). Given the classical description of any two ad- 
missible quantum channels Qo,Qi, the promised diamond norm problem PDN(Qo/ Q\> a , b) is 
asked to distinguish between the two cases , namely whether 1 1 Qo — Q\ \ | <> > a or | j Qo — Q\ \ | <> < b 
where a — b G Q(l /poly-log). 

As mentioned in the introduction part, Theorem Q] explicitly makes one conversion from dia- 
mond norms to equilibrium values that perverse the promised gap. It is easy to see that Theorem[T] 
works for any admissible channels Qo, Q\. Furthermore, this conversion can be done efficiently as 
long as the classical descriptions of Qo, Q\ are given. Thus, by combing the results in TheoremlH 
one can easily solve the promised version of diamond norm problems. Note that the input size 
changes to be the size of the matrix representing the channel now. 

Proposition 6. There is a NC algorithm that solves PDN(Qo, Q\, a, b) where a 2 — (4b — b 2 ) G 0(1 / poly-log) . 

Proof. This is a direct consequence when one combines the result of Theorem Q] and Theorem 0] 
Given the promise, by Theorem[TJ one can efficiently compute the equilibrium value A(H) whose 
value is either no more than v/4 2 ~ g2 or no less than 2 ^-. Hence if the difference a 2 — (4b — b 2 ) G 
Q(l/polylog), one can use the algorithm in Fig|2]to calculate A(E) efficiently in parallel to suffi- 
cient precision in order to distinguish between those two cases. The NC algorithm follows directly 
when one composes the circuits for each step. □ 

The only undesired thing of this algorithm is we can only solve the problem when the con- 
dition a 2 — (4b — b 2 ) G Q(l /poly log). This constraint makes it impossible for our algorithm to 
work for the whole range < b < a < 2. Evidences (implicitly in BRW05L e.g. Theorem 4.3) also 
demonstrate that simply repeating the channels for many times, like considering the alternative 
channels (Qo - Qi) m or Qf N - Qf N for some N, doesn't work either. 

Fortunately, there is one conceptually easy but technically detoured approach to amplify the 
gap in general. Particularly, we will make use of some known properties of the quantum inter- 
active proof systems and abuse them for a different purpose. Intuitively, we treat any quantum 
interactive proof protocol (assume the input is fixed) as a promised problem where the acceptance 
probability is either at least the completeness probability or no more than the soundness probabil- 
ity. Then we will convert the promised diamond norm problem into such a promised problem of 
one specific quantum interactive proof protocol. The crucial observation is the parallelization, am- 
plification lemmas together with the complete problem results about quantum interactive proof 
systems can be exploited to amplify arbitrary gap of any general promised diamond norm prob- 
lem and convert it to a new diamond norm problem^] that can be solved by Proposition|6j 

Let us demonstrate this approach with full detail. More importantly, we will show such con- 
version can also be computed efficiently. The latter one is due to the fact theparallelization and 
amplification procedures are constructed explicitly in [KWOOL Any protocol Iu with completeness 
a and soundness b will be denoted by P(a, b). Recall that any such protocol F(a, b) is treated as a 
promise problem where the acceptance probability is either at least a or at most b. The promised 
diamond norm problem PDN(Qo/ Q\,&,V) can thus be converted to the following protocol di- 
rectly. 

Definition 3. The protocol V [Qo, Qi]' 

6 Actually, it suffices to convert the original problem to close images problem. However, for the simplicity of descrip- 
tion, we choose QCD problem instead. 

7 Precisely, a single prover quantum interactive proof protocol. 
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1. If a 2 — (4b — b 2 ) G Q(l/polylog), use the algorithm in Proposition [6] to solve it directly. Oth- 
erwise, continue to next step. 

2. Convert the original problem to P [Qo, Qi] (1/2 + fl/4, 1/2 + b/4). 

3. According to the parallelization lemma and amplification lemma IKW00I , convert the proto- 
col P [Qo, Qi] (1/2 + a/4,1/2 + b/4) to P^(l, 1/2). According to the construction implicit in 
the proof of the QlP-Completeness of the problem QCDis,o.i, convert the protocol P^(l, 1/2) 
to a new promised diamond norm problem PDN(Q , Q[, 1.9,0.1). 

4. Use the algorithm in Proposition[6]to solve the new problem PDN( Q' , Q\, 1.9, 0.1) and return 
the answer. 



Figure 3: Algorithm for the PDN(Qo> Qi, a , b) problem. 

1. The verifier receives some quantum state p from the prover. 

2. The verifier selects {0, 1} uniformly and applies Qi to p and sends the result to the prover. 

3. The verifier receives some / from the prover, accepts if i = j and rejects otherwise. 

The protocol is almost identical to the Protocol 3.2 in IIRW05I . The only difference is the veri- 
fier needs to perform arbitrary admissible quantum channel Q,. It is not possible in general when 
the verifier's power is polynomial time bounded since arbitrary Q, might need huge number of 
gates to simulate. However, it won't be an issue for us since we treat such protocol as a promised 
problem without its original meaning and Qq, Q\'s description is already given. It follows im- 
mediately from IIRW05I that the protocol P [Qo/ Qi] has completeness 1/2 + a/4 and soundness 
1/2 + b/4 given the promise that either || Qq — Q\ || > a or || Qq — Q\ || < b. 

The parallelization and amplification lemmas IKW00| can then be reinterpreted as a way to 
convert any protocol P(c, d) to some protocol P'(c', d') with desired d and d! . This conversion can 
be efficiently computed when the gap between c and d is appropriate. For the promised diamond 
norm problem, we starts with some protocol P<>(fl, b) where the gap between a and b is at least 
inverse poly-logarithm and converts it to some protocol P^ (1,1/2). Moreover, such conversion 
can be done in NC. Secondly, because QCD\ 9 0.1 is QlP-Complete problem, we can convert the 
protocol P^(l, 1/2) again to a new promised diamond norm problem PDN(Qq, Q[, 1.9, 0.1) where 
such conversion is implicitly inside the proof of QlP-Completeness of the QCD\ 9 0.1 1IRW05H and 
the new channels Qq, Q\ can be computed in NC as well. Finally, we will invoke the algorithm in 
Proposition[6]to solve the new problem. To sum up, 

Proposition 7. There is a NC algorithm shown in Fig.\3\that solves general PDN( Qq, Qi, a, b) problems. 

As a standard technique, an algorithm for the promised version of problems can be used as a 
subroutine to solve the general problems via binary search when the range of the possible results 
is bounded. In our case, the diamond norm between any two admissible quantum channels is 
bounded between and 2. Hence, by recursively calling the subroutine in Fig. [3] O (poly-log) 
times, one can compute the diamond norm with Q(l /poly-log) precision. 
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Corollary 2. Given the classical description of any two admissible quantum channels Qq, Q\, the 
diamond norm of their difference \\Qo — Qi ||o can be approximated in NC with inverse poly- 
logarithm precision. 



6 Summary 

In this paper, we provide an alterative proof for QIP=PSPACE which starts from one QlP-Complete 
problem that computes the diamond norm between two quantum admissible channels. The key 
observation here is to convert the computation of the diamond norm to the computation of an 
equilibrium value. The later problem turns out to be a more structured problem and has a good 
solution in NC(poly) and thus in PSPACE. Besides reducing from the QCD problem, we could 
also reduce from the very first QlP-Complete problem close images MKWOOI or the protocol to sim- 
ulate QIP with competing provers in [GW05J. Both reductions will lead to the similar equilibrium 
values to the one in this paper. The technique of computing the equilibrium values in this paper 
can then be applied directly and lead to the same result. 

The multiplicative weights update method in our proof to solve the equilibrium value prob- 
lem can be generalized to solve a class of such equilibrium value problems. Particularly, for any 
density operator set D [X) and another convex compact set T, the following general equilibrium 
value problem 

A(O) = min max (c, ^>(p)) = max min (<&*(cr),p) 

can be solved efficiently in NC by the same algorithm in our paper if a good approximation algo- 
rithm to compute max^gp (^/^(i *)) given p* is available inNC and max^grj,^) max^r | (^/^(p)) I 
is bounded by some poly-logarithm function. 

One big open problem is to investigate to what extent the technique in this paper can be used 
to solve general equilibrium value problems. As mentioned in [JJUW09| , it is impossible to solve 
arbitrary SDPs in parallel unless NC=P. It might be the same case for the general equilibrium 
value problems. A recent effort HWulOl IGW10M made some progress on the general form of the 
equilibrium value that can be solved by similar techniques and the connection between the equi- 
librium value problems and semidefinite programming problems. Moreover, since the first re- 
lease of this paper, some deeper knowledge of variants of quantum interactive proof system is 
obtained. Particularly, the main open problem in | JJUW09| , namely whether QRG(2)=PSPACE, 



is resolved MGW10II with positive answer. The class QRG(2) contains all problems which can be 
recognized by two-turn (i.e, one-round) quantum refereed games. The classical analogue of this 
class is known to coincide with PSPACE [FK97H. 

Another open problem is how to extend the connection we build here between the computa- 
tion of diamond norms and the computation of equilibrium values. One might hope to obtain 
efficient parallel algorithm for calculating diamond norms of any quantum channel. 
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A Proof of Theorem |3] 



Assume < < 1 for all t, after T rounds, the algorithm in Fig [TJ guarantees that, for any 

p* G D(X), we have 



(l-e)£(pM MM)<(p*,X> 



r(0 



InN 



f=l \ i=\ I 

Proof. It is easy to see that all WW G Pos (X). Observe that, for t = 1, • • • , T, 



Tr(W (m) ) =Tr 
= Tr 



exp ( -e £ m(t) 
W w exp (-eM (f) 



< Tr 
< Tr 



exp y-e £ M (t) J exp ( -eM w 
W®(l x -e!M®) 



(W®,l x - e'M«) = Tr [w«l (1 - e' (p« M«)) 
WW] exp(-£ / ^W,MW)) 



< Tr 



The inequality in the first line is due to Golden-Thompson inequality [Bha97|. The second in- 
equality is due to Lemma [8] where e' = 1 — e~ £ . The third line is obtained by substituting pW = 
WW/TrWW. The final inequality is obtained by noticing that 1 — e'x < e~ ex ,"ix G [0,1] and 

'p(t),M®) G [0,1]. 



If we repeat the process above, by induction as well as the fact IAA 1 ) = \%, we have: 



Tr 



W( T+1 )j < Nexp(- £ / £ (p^,M^ 



T=l 



where N = dim (X) as defined above. On the other hand, we have 



Tr 



Tr 



>exp(-eA N (£MW)) 



T=l 



T=l 

The last inequality holds because Tr(e A ) = ^ e A <( A ) > e A i( A ). Thus, we conclude that 

T T 

exp(-eA N (£ MW)) < Nexp(-e' £ (pW,M< f) )) 



T=l 



T=l 



Since is the minimum eigenvalue, for any density operator p* G D ( X ), we have Ajv (Xt=i ^ ) — 
<^j0*, Et=i V Take the logarithms of the both side and simplify as well as notice the fact s! ob- 
tained by Lemma[8]has the property that e' > e(l — e), then we have: 



(l-e)f)(pW # MW)</^EM« 



InN 



for any density operator G D (A") as required. 



□ 
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For the sake of completeness, we prove the lemma we used in the proof below. 

Lemma 8. For any < e < 1/2, let e' = 1 — e~ £ . Then we have the following matrix inequality, for any 
< M < 1, 

exp(-eM) < 1-e'M 

Moreover, we have e' > e(l — e). 

Proof. It is easy to verify that for any e G [0, 1/2], we have 

f{x) = exp(— ex) < g(x) = 1 — e'x if iG [0,1] 

Since < M < 1, let M = UDU f be the diagonlization of M. Then,/(M) - g(M) = Lf(/(D) - 
g(D) ) IT^. Since D is a diagonal matrix of which every diagonal entry contains one eigenvalue of M 
and thus is in [0, 1], then /(D) — g(D) is a diagonal matrix with non-positive diagonal due to the 
inequality above. Thus, /(M) — g{M) is a negative semidefinite matrix and hence / (M) < g(M), 
namely, exp(— eM) < 1 — e'M. It is also easy to verify that e' > e(l — e). □ 

B Comments on precision issues 

The analysis made in the main part of this paper has assumed that all computations performed by 
the algorithm are exact. However, in order to implement our algorithm, some step of the compu- 
tations must be approximate. Particularly, the computation of the positive eigenspace projections 
and the matrix exponentials will need to be approximate. An elaborated analysis on these issues 
can be found in |JJUW09| |JW09 [ . We will basically follow that type of analysis and provide a brief 
sketch of the analysis to the specific problem in our paper. 

First, it must be made clear which part of the algorithm can be made exact and which part 
must be made approximate. We will use the same convention of storing complex numbers as the 
one in |JJUW09| . Once the input x is given and stored in memory, all elementary matrix opera- 
tions (in this case: addition, multiplication, and computation of the trace or partial trace) can be 
implemented exactly in NC MGat93L However, the matrix exponentials and positive eigenspace 
projection cannot be exact since these operations will generate irrational numbers and the preci- 
sion must be truncated at somewhere. Fortunately, Watrous et al. [JJUW09| provided a way to 
approximate those two operations to high precision in NC. Precisely, 

Fact 2. Given an n x n matrix M (whose operator norm bounded by k) and a positive rational 
number rj, the computation of n x n matrix X such that || exp(M) — X \\ < n can be done in NC. 

Fact 3. Given an n x n Hermitian matrix H and a positive rational number n, the computation of 
an n x n positive semidefinite matrix A < 1 such that || A — A|| < n for A being the projection 
operator onto the positive eigenspace of H can be done in NC. 

Before we move on to the analysis of the precision issue, it helps to introduce the following 
convention. We will represent the actual matrices generated during the algorithm by placing a 
tilde over the variables that represent the idealized values. As we discussed above, there are 
mainly two types of operations where the accuracy will be lost. Further investigation tells us 
that the matrix exponentials are always necessary to the multiplicative weight update method 
while the positive eigenspace projections are special for our application. For the generality of the 
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analysis, we will first discuss what the general form of Theorem[3]is when the computation is only 
approximate. 

Consider the scheme in Fig [TJ and keep the notation convention in mind. The will be the 
actual generated density operator for round t and w( f+1 ) = exp(—e J^=i M w ). The latter one is 
exact simply because WW is only a notation and not stored in the memory at all. Fact[2]implies that 
||p(f) _ w( f )/TrW( f )|| <<S!/Nfor every t where 5\ is some constant for our purpose. The situation 
for is tricky in the sense that there is no idealized value for to have in the general scheme. 
By going through the proof of Theorem [3] again, we can easily obtain the following fact. 

Fact 4. If the computation can only be performed approximately, the inequality in Theorem [3] 
becomes 

(1 - e) t (p« MW ) < L £ m \ + l lf + \ T5l 

Now consider the concrete in Fig|2j By making use of the fact above we can repeat almost 
all the steps in the proof of TheoremHl The only change is to replace the Equation Q by 




By Fact|H we have \\Yl^ — Yl^> \\ < 5\/N where the Tv?) is the projection onto the positive 
eigenspace of S(p w ). Please note that 11^ here is not its idealized value when everything is exact 
but rather the exact value given the approximate p^\ The rest part of the proof follows similarly. 
Finally, we will get 

T 

A(H) - St < A = i £ (pM, E* (flW)) < A(E) +S + S 1 

Since our target is to distinguish between two promises with constant gap, we can choose 8\ 
to be any sufficiently small constant. In this way, the precision issues are handled. 
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